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SUMMARY 

A steady, three-dimensional average-passage equation system Is derived for 
use In simulating multistage turbomachinery flows. These equations describe a 
steady, viscous flow that Is periodic from blade passage to blade passage. 

From this system of equations, various reduced forms can be derived for use In 
simulating the three-dimensional flow field within multistage machinery. It Is 
suggested that a properly scaled form of the average-passage equation system 
would provide an Improved mathematical model for simulating the flow In multi- 
stage machines at design and. In particular, at off-design conditions. 


INTRODUCTION 

The flow field In multistage compressors and turbines Is extremely com- 
plex. This flow field Is highly unsteady, with time scales ranging from a 
fraction of shaft speed to several times that of the highest blade passing 
frequenrv. The length scales of the flow arc also diverse. They range from 
the circumference of the machine to a fraction of a blade chord. The vast 
range of time and length scales makes direct numerical simulation practically 
Impossible. Confronted with analyzing complex, multiple-scale problems, aero- 
dynamlclsts have traditionally described their physics In terms of an appro- 
priately averaged set of equations. A prime example of this methodology Is the 
Reynolds-averaged modeling of turbulent flows. For turbomachinery the litera- 
ture seems to support a description based on an average blade passage. With 
respect to a given blade row the flow within the average passage Is steady and 
spatially periodic from blade passage to blade passage. 

In the present work a rigorous attempt Is made to derive the equations 
governing this flow field. To the author's knowledge, this three-dimensional 
equation system has not been defined. It can be reduced to the familiar com- 
bined through-flow (meridional), blade-to-blade (cascade) equation set (ref. 1) 
by Introducing a number of assumptions, all of which appear to be justified for 
highly efficient machines operating near design. Further, they can be reduced 
to the through-flow equations derived by Sehra (ref. 2) by assuming the flow to 
be Invlscld, by assuming that heat Is added reversibly, and by employing a tan- 
gential averaging operator. More Inportantly, the derived equation set pro- 
vides a reference against which multistage analyses .can be judged for their 
completeness. In particular, given the advances being made In the field of 
computational fluid dynamics. It might be time to determine If any of the 
assumptions that restrict the validity of today's multistage analyses to the 
neighborhood of the design point can be removed. Of these assumptions those 
that are consistent with a quasl-three-dlmenslonal flow structure need to be 
carefully reevaluated. If these assumptions can be removed without yielding 
an equation system that Is essentially unsolvable on second-generation class VI 


computers, the resulting equation system should enhance our ability to simulate 
the flow In multistage turbomachines. 

Three averaging operators are developed for application to the Navler- 
Stokes equations In the next section of this report. The first averaging oper- 
ator, referred to as 'ensemble averaging," Is Introduced to eliminate the need 
to resolve, In detail, the structure of turbulent flows. The application of 
this operator to the Navler-Stokes equations yields the Reynolds-averaged form 
of these equations. This equation set describes an unsteady deterministic flow 
field. The global effects of turbulence on this flow field are accounted for 
by means of a Reynolds stress and a scalar flux tensor. 

The second operator to be developed Is a time-averaging operator. This 
operator Is used to average the deterministic unsteady equations In time at 
every point In space. This operation removes unsteady time scales that are of 
the order of the period of shaft rotation or less. The global effect of the 
coherent (as opposed to random) blade-to-blade unsteady flow structure mani- 
fests Itself In the resulting equation system as body forces, energy sources, 
momentum, and energy temporal mixing correlations. The effect of these terms 
on the axlsymmetrlc flow field generated by an Isolated rotor was Investigated 
by Sehra (ref. 2). For a single-stage machine subjected to an axlsymmetrlc In- 
let and exit condition, the resulting equations describe the three-dimensional 
average-passage flow field associated with either the rotor or stator blade 
row. For multistage machines the equation set governs. In general, a flow 
field that varies from passage to passage around a given blade row. 

For the average-passage equations to be extracted from this equation 
system, they must be averaged on a passage-to-passage basis. The third averag- 
ing operator that Is developed accomplishes this task. This operator averages 
out the details of the passage-to-passage variation In the flow field. How- 
ever, the global effect of this variation on the average-passage flow field Is 
not eliminated. Its existence Is accounted for through body forces, energy 
sources, momentum, and energy spatial mixing correlations. 

The average-passage flow equations are derived In the third section of 
this report. These equations contain correlations that arise from turbulence, 
unsteady coherent flow, and passage-to-passage flow variations. These correla- 
tions are generic to this equation system, as the Reynolds stress tensor Is 
generic to the Reynolds-averaged Navler-Stokes equations. They serve as the 
mechanism by which the passage-average flow field Is energized In a multistage 
environment. The need to account for this mixing process In turbomachinery 
flow analyses was recently documented In reference 3. 

Two reduced forms of the average-passage flow equations are developed In 
the final section of this report. In deriving the first set of equations It Is 
assumed that the body forces and energy sources generated by the presence of 
neighboring blade rows can be evaluated In a through-flow blade-to-blade anal- 
ysis (ref. 4). It Is further assumed that the mixing correlations can be mol- 
eled In both the through-flow blade-to-blade analysis and the average-passage 
system according to the analysis presented In reference 3. The resulting set 
of equations describes a three-dimensional flow field whose axlsymmetrlc com- 
ponent Is equal to that predicted by the through-flow blade-to-blade model. In 
the second approach an average-passage equation system Is developed for each 
blade row In the machine. The resulting system of equations Is coupled through 
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common expressions for the mixing correlations, the body forces, and the energy 
sources. Thus the equations must be solved simultaneously. The closure prob- 
lem associated with this formulation requires the development of models that 
globally describe the mixing produced by the coherent spatial and temporal non- 
uniform flow. In addition to models for the Reynolds stress and associated tur- 
bulent correlations. 


AVERAGING PROCEDURE 
Ensemble-Averaging Procedure 

To derive the average-passage equation system for turbomachinery flows, 
three averaging operators must be applied consecutively to the Navler-Stokes 
equations. The first operator, often referred to In the literature on turbu- 
lence as "ensemble averaging" (ref. 5), serves to decouple the unsteady deter- 
ministic flow features from the random fluctuations. This averaging procedure 
yields the familiar Reynolds-averaged Navler-Stokes equation. Its mathematical 
form Is defined as 
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(1) 


where f^ denotes the 1^ realization of the variable f and N denotes the 
total number of observations of this variable. That Is, If N experiments are 
conducted, each experiment produces a realization (l.e., an observation) of a 
physical variable at a given point In spa^e at a given Instant In time. The 
average value of this variable Is estimated by evaluating the sum of the reali- 
zation divided by the number of experiments. It Is assumed that the quotient 
becomes Independent of the number of experiments as the number of experiments 
becomes large. 

For the analysis of a compressible turbulent flow It has proven useful to 
define a density-weighted ensemble average (ref. 6). The form of this averag- 
ing operator Is 


? 


11m 

N-kd 




Vi 


1=1 


( 2 ) 


where p Is the ensemble average density end p^ Is the 1*h realization of 
the density field. The velocity field as well as the fluid total properties Is 
generally defined In terms of this weights average. Finally note that this 
averaging operator commutes with differentiation; that Is, 


af af 
dx m ax 


(3) 
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where x Is an appropriate Independent variable. This property Is often taken 
for granted In deriving the Reynolds-averaged Navler-Stokes equations. 


Time-Averaging Procedure 


The second averaging operator yields the time-averaged flow equations. 

This operator Is applied throughout the flow domain, Including the region with- 
in the rotating blade rows. Several researchers have used a time-averaging 
operator to develop the steady flow equations for the flow field downstream of 
a rotating blade row (ref. 7). In the current analysis the range of applica- 
bility of this operator Is extended to Include the flow region within blade 
rows (l.e., stationary and rotating). The traditional definition of this aver- 
aging operator Is given by the Integral 
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where T Is the period of one rotor revolution, t. Is the reference starting 
time, and F the deterministic component of f. To extend this averaging pro- 
cedure to the flow region within the rotating blade rows, a "gate function" 
must first be defined for each rotating blade row. This function Is represent- 
ed by a finite number of step functions. Its value Is 1 for any point within 
the flow domain and zero In the regions occupied by the blades. The mathemati- 
cal form of this function for a blade row with N blades Is 
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where U Is the unit step function, ft Is the rotational speed of the shaft, 

Kp Is the number of shaft revolutions over which H Is defined, and t re f Is 
tne reference starting time. The remaining variables In this equation are de- 
fined In figure 1. The length of time during one revolution of the shaft for 
which this function Is nonzero Is 
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where (02 - ®i) defines the local blade thickness. The ratio of this length 
of time to the period of revolution (T = 2 p/Q) Is a measure of the blockage 
due to the blading (ref. 1). We define for future us»> a weighting function \ R 
that Is equal to 
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(l.e., 1 minus the blockage) within a rotor blade passage and 1 outside a 
passage. 

A time-averaging operator that Is applicable to all regions of space can 
be defined as 
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H(t)f(t)dt 
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Outside a rotor passage this operator reduces to equation (4), since Xp and 
H become equal to 1. Under these circumstances It can be shown that the time- 
averaging operator commutes with differentiation. However, within a rotor pas- 
sage this property Is lost. 


A formula that specifies the rules governing the order of spatial differ- 
entiation and time averaging at every point in space Is needed In the next 
section. This formuli will be developed for the r partial derivative a/a r . 
The formulas for the partial derivatives with respect to either o or Z can 

be developed In a similar fashion. Begin by directly substituting df/ar Into 
equation (8) to get 
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This expression can be rewritten as 
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since the limits of Integration are Independent of r. The two Identical In- 
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tegrals In equation (10) are equal to f. The remaining Integral can be eval- 
uated by noting that 
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where 4 Is Dirac delta function (ref. 8). Introducing this expression Into 
the last Integral In equation (10) yields 
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where t Is an Integer. Equation (10) can be rewritten as 
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This equation Is the formula for Interchanging time averaging and differentia- 
tion with respect to r. It Is valid at every point In space. A similar re- 
sult can be derived for Interchanging differentiation w*th respect to either 
e or Z and time averaging. A formula for Interchanging time averaging and 
the temporal derivative will also be needed In the next section. To derive 

this formula, substitute af/at Into equation (8) to get 
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This expression can be Integrated by parts to yield 
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According to the Leibnitz rule (ref. 9), the quantity within brackets can be 
rewritten as 
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The Integral In equation (15) can be evaluated once the expression 
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Is substituted for aH/at. Thtf final result Is 
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Note that outside a rotating blade passage the expression within brackets In 
equations (13) and (18) vanishes and that x becon.es equal to 1. Hence time 
averaging commutes with either temporal or spatial differentiation outside a 
rotating passage. 

For a compressible fluid It will prove useful to define the velocity and 
total flow properties In terms of their density-weighted time average. This 
average Is defined as 
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which Is analogous to the definition of the density-weighted ensemble average. 
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Passage-to-Passage Averaging Procedure 


For a single-stage machine or a repeating-stage multistage machine sub- 
jected to axl symmetric Inflow and outflow, the flow field described by the 
steady-flow equations will be periodic from blade passage to blade passage. 

Thus these equations, by definition, govern the average-passage flow field. 
However, for a multistage machine In which the number of blades varies from row 
to row, the steady flow will not, In general, be Identical In each passage of 
a given row. For the average-passage flow equations for a general configura- 
tion to be developed, a third averaging operator must be Introduced. This 
operator will be developed by means of the theory associated with Fourier anal- 
ysis. This development Is begun by defining a gate function for each stator 
blade row. (Observe that the analysis being developed Is for the absolute 
frame of reference.) This function, like the preceding one (l.e., eq. (5)) Is 
defined In terms of a finite number of unit step functions. For a stator blade 
(rotor for the relative frame of reference) row with L blades this function 
Is equal to 
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The variables a, 63 , and 64 are defined In figure 2. The ratio of the 
angular distance for which G Is nonzero to the angular distance around the 
machine (l.e., ?«) Is a measure of the blockage attributed to the blade row. 
This ratio Is equal to 
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The gate function G Is used In the construction of a Fourier series 
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representation In e of the flow variable f. This representation Is given 
by the equation 
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where the Fourier coefficients A n are defined as 
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The passage-to-passage average of f Is defined with respect to a referenced 
stator (rotor) blade row. For a blade row with M blades this average Is de- 
fined as 
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where A k || Is the kM Fourier coefficient and the weighting function x§ Is 
set equal to 1 outside a stator passage and for all stators having multiples of 
N blades. For the remaining stators x§ Is given by equation (21). By con- 
struction, f Is periodic over a distance of 8 » 2*/M. Introducing the Inte- 
gral representation for A kM (l.e., eq. (23) Into eq. (24)) yields 
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Interchanging the order of summation and Integration In this equation and not- 
ing that the Infinite series can be expressed as (ref. 8) 
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The order of Integration and summation Is reversed to simplify this result. 

The Integral that appears under the summation sign can then be evaluated. The 
final form of the expression for f becomes 
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Equation (28) states that the passage-to-passage average of a variable Is equal 
to the arithmetic mean of the value of the variable at the corresponding rela- 
tive position within a blade row times the blockage factor 1/XM. Thus the 
passage-to-passage average of a variable will In general be different for each 
blade row. It will, however, be periodic In e, the period being set jy the 
geometry of the referenced blade row. 

In the next section the passage-to-passage averaging operator Is used to 
transform the steady-flow equations Into their passage-to-passage averaged 
form. Before carrying out this exercise. It Is necessary to establish the 
rules governing Interchanging passage-to-passage averaging and spatial differ- 
entiation. This rule Is developed here for the r partial derivative. The 
rule for the remaining spatial derivatives (l.e., 8 and z) can be developed 
accordingly. Equation (25) states that the passage-to-passage average of 

af/ar with respect to a blade row with M blades Is 
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Integrating this expression by parts yields 
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The expression within brackets In equation (30) Is equal to ? (refer to 
eq. (25)). To evaluate the Integral that appears In the last term In 
equation (30), Insert for dG/dr the expression 



Carrying out the Integration yields 
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This result can be rewritten as 
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The first terra In this expression represents the axlsyrametrlc contribution; the 
remaining two terras represent the nonaxl symmetric contribution. Introducing 
this result along with equation (25) Into equation (30) yields 



This equation represents the formula for Interchanging the order of passage- 
to-passage averaging and partial differentiation with respect to r. Similar 
expressions can be developed for the 6 and z partial derivatives. The 
partial derivative with respect to t-j can be shown to commute with the 
passage-to-passage averaging operator. 

For a compressible fluid It Is advantageous to define the velocity field 
and total flow variables In terms of a density-weighted average. The density- 
weighted form of the passage-to-passage averaging operator Is defined as 
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This definition Is In accordance with equation (19). 

This completes the definition of all of the averaging operators needed to 
construct the average-passage equation system for a multistage machine. These 
definitions, along with the supporting analysis, are used In the next section. 


PASSAGE -AVERAGEO FLOW EQUATIONS 

The equations describing the average-passage flow field are derived by 
sequentially applying the three averaging operators defined In the previous 
section to the Navler-Stokes equations. The continuity and tangential (l.e., 

6) momentum equations associated with this flow field are derived here. The 
remaining two momentum equations (l.e., r and z direction) and the energy 
equation are given In the appendix. 

To begin, the velocity field (defined In the absolute frame of reference) 
Is decomposed Into Its density-weighted ensemble average component plus a ran- 
dom variable. The form of this decomposition for the radial velocity component 
V r is 
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where the superscript denotes the random component. A similar expression ex- 
ists for the tangential velocity component V 8 and the axial velocity com- 
ponent V z . The continuity equation expressed In terms of the cylindrical 
coordinates r r e, and z Is 
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Introducing the suggested velocity decomposition Into equation (37) and multi- 
plying the result by the ensemble-averaging operator produces 
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slrce the density-weighted ensemble average of the random velocity field Is 
zero. 


Tr. next step entails decomposing the deterministic velocity field Into 
Its density-weighted, time-averaged component plus an unsteady component. For 
the radial component of velocity, this composition Is of the form 


V r - V r ♦ \ (39) 

where the superscript * denotes the unsteady component. A similar decomposl 
tlon will be Introduced for the two remaining velocity components. Introducing 
equation (39) along with the corresponding expressions for the ether velocity 
components Into equation (38) yields, upon time averaging, the result 
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where Xr Is the weighting function associated with the blockage of the rotat- 
ing blade rows (eq. (7)). In deriving equation (40) the rules for reversing 
the order of differentiation and time averaging as derived In the preceding 
section were used. 


The final step entails decomposing the steady velocity field Into a 
density-weighted, passage-to-passage averaged component and an aperiodic com- 
ponent. Decomposing the steady radial velocity component In this fashion 
yields 
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where the superscript ~ denotes the aperiodic component. Introducing this 
decomposition for each velocity component Into equation (40) and evaluating Its 
passage-to-passage average with respect to the J th stator yield 
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where \\ Is a weighting function associated with the blockage of the neigh- 
boring blade rows, ‘ibis function Is equal to 1 outside a blade passage. With- 
in a rotor or stator passage It Is set equal to Xr or x$, respectively. 

This equation follows directly from the material presented In the preceding 
section. Equation (42) represents the continuity equation for the average- 
passage flow field. Its mathematical structure Is very similar to the conti- 
nuity equation for the Navler-Stokes system. 

The tangential momentum equation In the absolute frame of reference Is 
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where T r e, t ee , t z @ are components of the viscous shear stress and p Is the 
static pressure. The derivation of the passage-to-passage average form of this 
equation is begun by Introducing the velocity decomposition, as defined by 
equation (36) for all three velocity components. Into equation (43). Evaluat- 
ing the ensemble average of the resulting expression yields the Reynolds- 
averaged form of the tangential momentum equation. This equation has the form 
(ref. 10). 
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; nders J ored ^rras represent three of the nine components of the 
Reyno^s stress tensor. The time-averaged form of this equation Is derived bv 
first decomposing each velocity component according to equation (39) The re* 
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The algebra that led to this result made use of the material Dresented in th* 
preceding section and the velocity product averaging rule 
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(Similar expressions can be derived for the products P V V and pV V .) The 

underscored terms In equation (45) represent the mixing stress attributed to 

the unsteady coherent velocity field. The body forces f[® R) and F (eR) are 

associated with the Invlscld and viscous forces acting on the rotating blades 
The Invlscld component F is equal to rotating Diaaes. 
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The corresponding expression for the viscous force 1: 
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Both of these forces are axlsymmetrlc and vanish outside a rotor passage. 


( 47 ) 


For an Isolated stage equation (47) Is the tangential momentum equation 
of the average-passage equation system. For this limiting case the aerodynam- 
ic loading generated by the rotor produces an axl symmetric body force and an 
axlsymmetrlc energy source, which are distributed over the axial extent of the 
rotor passage. The unsteady coherent velocity field generated by the rotor 
produces a stress that acts throughout the flow domain. 

As was noteu previously the time-averaged flow field within a nultl stage 
machine Is not. In general, periodic from blade passage to blade passage. The 
steady tangential momentum equation, like the steady continuity equation, can 
be transformed Into Its average-passage form by using the passage-to-passage 
averaging operator. This transformation Is Initiated by decomposing the veloc- 
ity field In equation (45) according to equation (41). The resulting equation 
Is then operated upon by equation (25). The reference blade row for this oper- 
ation Is the stator. The resulting expression for the passage-to-passage 

averaged tangential momentum equation Is of the form 
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Similar expressions can be derived for the radial and axial momentum equations. 
They are given In the appendix, along with the passage-to-passage averaged 

energy equation. The variables F and Fy 6 ^ In equation (48) represent 

the Invlscld and viscous body forces generated by the neighboring stator blade 
rows. They are equal to 
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The expressions In parentheses In equations (49) and (50) 
metric component of the body force; the remaining terms represent the nonaxl- 
symmetrlc component. These forces vanish outside the stator passages. The 
underscored terms In equation (48) represent the mixing stress generated by the 
aperiodic component of the steady velocity field. Thus the total mixing stress 
associated with the average-passage equation system Is composed of a stress due 
to the steady, aperiodic velocity field, a stress generated by the coherent. 


¥) 


(50) 

represent the axlsym- 
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unsteady velocity field, and the Reynolds stress. The total stress can be ex- 
pressed as 


A X 

R U ♦ pMj (51) 

where subscripts 1 and J are assigned the values 1, 2, and 3, which corre- 
spond to the radial, tangential, and axial coordinates, respectively. This 
second-order tensor Is generic to the average-passage equation system. Its 
evaluation, along with the body forces and energy sources, constitutes the 
closure problem for this equation system. Two possible approaches for solving 
this problem are presented In the next section. 

This completes the theoretical development of the average-passage equation 
system. This equation system was derived without Introducing any additional 
assumptions other than those encompassed In the Navler-Stokes equations them- 
selves. The equation system was constructed to be spatially periodic, with a 
period equal to the pitch of a reference blade row. Hence, In general, there 
Is one such equation system for each blade row of a multistage machine. 


CLOSURE MODELING 

The average-passage equation system, like the Reynolds-averaged Navler- 
Stokes equations, does not contain sufficient Information to determine Its 
solution. To close this system of equations, one must construct a mathematical 
model for the mixing stress tensor R^ 'eq. (51)), the body forces, and the 
energy sources. Two approaches to this problem are outlined. 

The first approach makes use of the capability that exists for simulating 
the axlsymmetrlc flow field In multistage machines (ref. 4). This analysis, 
often referred to as "through -flow analysis," has been tuned over the years to 
accurately simulate the flow In turbomachinery at design and near-design oper- 
ating conditions. This form of analysis Includes models for the body forces 
and the energy sources. These same models could also be used In the present 
analysis. 

Recently a model has been proposed by Adkins and Smith (ref. 3) by which 
the effects of mixing can be Included In through-flow analyses. They suggest 
that the angular momentum, total temperature, and total pressure field be mixed 
according to a simple two-dimensional diffusion equation. The diffusion co- 
efficient associated with this equation Is developed from Invlscld, secondary- 
flow theory that has been empirically corrected for the effects of viscosity 
and Interactions with downstream bladings. Their work could serve as a basis 
for developing a model for the mixing stress as defined In the present analy- 
sis, since It Is this stress field that produces the mixing they attempted to 
model. Away from solid boundaries the diffusion model Is analogous to 
Bousslnesq's eddy-viscosity model. Thus It appears reasonable to assume that 
R^j can be approximated by an eddy viscosity model In the Interior flow 
region, the eddy viscosity being directly proportional to the diffusion coeffi- 
cient of reference 3. Near solid boundaries this model must be modified to 
acknowledge the presence of the boundary. This modification could resemble 
that used for the Inner region of a turbulent boundary layer. The resulting 
mixing stress model would thus be a composite. A similar approach, based on 

17 


Bousslnesq's eddy-conductivity model, could be used to develop expressions for 
the correlations Involving the product of velocity and enthalpy that appear In 
the energy equation (A3). This approach, along with the suggested models for 
the mixing stress, the body forces, and the energy source, Is sufficient to 
close the average-passage equation system. 


An alternative approach to the closure problem Is to address separately 
the contribution to the mixing stress, the body forces, and the energy sources 
made by the coherent velocity from that made by the turbulent fluctuations. 

The contribution of the coherent velocity field to the mixing correlations can 
be decomposed Into that attributed to an Incident gust and that attributed to 
a radiated velocity field. The radiated field results from the Interaction of 
the gust with blading. The Incident-gust velocity field Is constructed by 
superimposing the passage-to-passage averaged velocity field of each blade row. 
Its mathematical form referenced to the n™ stator Is 


J 


max 


V (G) - 
V 1 " 



J=1 



(r,e,z) - V L (r,z) 


max 


E 


k=l 


v (k) (r,e - at.z) 


V^(r,z)| 


!= 1,2,3 (52) 


where V, denotes the mass-averaged axlsymmetrlc velocity field, 0 the 

1 ma X 

number of stationary blade rows, and K the number of rotating blade rows. 

max 

Added to this velocity field Is a radiated field that, because of the boundary 

conditions applied at solid surfaces, must be equal and opposite to ' at 

these surfaces. An evaluation by numerical simulation of this field without 
Introducing any simplifying assumptions would require a direct simulation of 
the unsteady deterministic flow field In a multistage machine. Such a simula- 
tion Is well beyond the limits of any of today's supercomputers. Including the 
Cray II. Therefore semlemplrlcal methods must be resorted to for modeling the 
radiated velocity field and Its contribution to the various correlations. For- 
tunately the problem of predicting this field has been addressed In a number 
of turbomachinery research programs over the years. In particular the problems 
of forced vibration and aerodynamic noise require a knowledge of the radiated 
velocity field before fluctuating airloads and noise can be predicted. A num- 
ber of models have been developed as a result of research In these problem 
areas. Because these models are In general only applicable to the Invlscld 
flow region, they must be extended to Include the viscous region surrounding 
the blading. This will be a formidable task. It will require detailed exper- 
imental measurements and highly refined numerical simulations to reveal the 
physics of the unsteady-flow process In this region. The data obtained from 
these studies will provide a basis for developing models for the correlations 
associated with the deterministic velocity field. Parallel with this effort, 
a research program Into the effects of a highly unsteady coherent free stream 
on the turbulent boundary-layer correlations (l.e., Reynolds stress, etc.) 
needs to be pursued. The unsteady coherent structure could have a significant 
effect on the turbulent correlations since the length and time scales of coher- 
ent structure can be comparable to the length and time scales of the large 


eddies In a turbulent boundary layer. Finally the Importance of the turbulent 
correlation relative to those generated by the deterministic velocity field 
will have to be assessed. 


CONCLUDING REMARKS 

This report presents a rigorous derivation of the average-passage equa- 
tion system. These equations have a closure problem that Is quite similar to 
that for the Reynolds-averaged Navler-Stokes equations. The task of scaling 
the terms In this equation set has not been undertaken. If, however, this task 
could be executed successfully and If the closure problem associated with the 
dominant correlations could be solved, the resulting scaled system could pro- 
vide a more fundamental basis for the analysis of turbomachinery flow than cur- 
rently exists. 


APPENDIX - AVERAGE-PASSAGE EQUATION SYSTEM 


This appendix summarizes the three-dlmenslona.1 average-passage equation 
system. The continuity equation and the tangential momentum equation were 
derived In the main body of this report. They will be restated to complete the 
equation set. The governing field equations averaged with respect to the 
stator are 

Continuity: 
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The blade forces F ( jjj S) , F* RS) , F ( j RR) , and F< RR) are defined as 
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Tangential momentum equation: 
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The blade forces 
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Axial momentum equation: 
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Energy equation 
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(All) 
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where the energy source E Is equal to 
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In equations (All) and (A12) e 0 Is the total Internal energy, H 0 the 
total enthalpy, K the thermal conductivity, and T the absolute temperature. 
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Figure 1. - Starting configuration of a rotor b>jde row. 





